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I give a brief review about the color glass condensate, which is the universal form of hadrons and 
nuclei at high energies. 



1 The Color Glass Condensate 

The main purpose of this talk is to con- 
vince you that high energy limit of QCD is 
the Color Glass Condensate. 1 More precisely, 
when the energy of hadronic scattering is very 
large, the participating hadrons or nuclei be- 
have not as usual states made of valence par- 
ticles, but as new form of matter called the 
color glass condensate (CGC). This name is 
after the following observations. First of all, 
it is made of gluons (" small- a; gluons" ) which 
have color and carry small fractions x <C 1 
of the total momentum. Next, these small- 
x gluons are created by slowly moving color 
sources (partons with larger x) which are dis- 
tributed randomly on the two dimensional 
disk (the Lorentz contracted hadron). This is 
like a glass whose constituents are disordered 
and appear to be frozen in short time scales. 
Lastly, the density of small- a: gluons becomes 
very large until it is saturated to some value. 
Typically the occupation number of gluons 
is of 0(l/a s ) > 1 at saturation, which is 
like a condensate of bosons. As the scatter- 
ing energy is increased, the hadrons undergo 
multiple production of small- a; gluons, and 
eventually become the CGC. 

Note that the claim above is as correct as 
the statement about the other limits of QCD: 
high temperature/density limit of QCD is 
the QGP/color superconductor, which now 
everyone believes true. In the same sense, 
if one goes to high energy limit in QCD, 
one will necessarily encounter the CGC. Note 
also that these three different limits allow for 
weak-coupling descriptions powered by so- 
phisticated resummation schemes. 



2 Gluon saturation and unitarity 

Let us explain how the CGC appears with 
increasing energy. When the gluon density is 
not so high, change of the gluon distribution 
with increasing energy (or decreasing x) is 
described by the BFKL evolution equation: 2 
d 

—N T (k)=a s K BFKL ®Af T (k), (1) 

OT 

where a s = a s N c /ir, N T (k) is the dipole- 
proton scattering amplitude (or unintegrated 
gluon distribution of the target proton), r = 
lnl/x is the rapidity, and <5 s -/ybfkl is the 
kernel representing the probability of split- 
ting of one dipole into two. This is essentially 
a linear differential equation, and its solution 
at asymptotically large energies shows expo- 
nential growth of N T (k). Namely, multiple 
production of gluons occurs endlessly. This, 
however, violates the unitarity bound for the 
cross section (or, N T < 1) and the BFKL 
equation must be modified so as not to vi- 
olate the unitarity. In fact, what is missing 
in the BFKL equation is the recombination 
process of two gluons into one, which cannot 
be ignored when the gluon density is high. 
Note that this process effectively reduces the 
speed of growth. Once this is included, the 
BFKL equation is replaced by the Balitsky- 
Kovchegov (BK) equation: 3 

jV r (fc) = a s K BPKL ®(Af T (k) - AA T 2 (fc)). (2) 

Compared to the BFKL solution, the solution 
to this equation is drastically changed due to 
the presence of the nonlinear term. Indeed, 
it exhibits saturation (unitarization) of the 
amplitude and a kind of universality, as I will 
explain below in a simple example. 
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2.1 Analogy with population growth 
In order to understand what happens in the 
BK equation, let us ignore the transverse dy- 
namics for the time being. This simplifica- 
tion allows us to find an interesting analogy 
with the problem of population growth. Long 
time ago, Malthus discussed that growth rate 
of population should be proportional to the 
population itself, and proposed a simple lin- 
ear equation for the population density N(t): 



-N(t) = aN(t). 



(3) 



Its solution N(t) = N Q e Qt shows exponen- 
tial growth known as the "population explo- 
sion." However, as the number of people in- 
creases, this equation fails to describe the 
actual growth. This is because various ef- 
fects such as lack of foods reduce the speed of 
growth. One can effectively represents such 
effects by replacing the growth constant a by 
a(l — N) which decreases with increasing N. 
This yields the famous logistic eguation first 
proposed by Vcrhulst: 

±N(t) = a(N(t)-N 2 (t)). (4) 

This nonlinear equation can be solved ana- 
lytically. The similarity of these equations to 
our problem is rather trivial: N and t corre- 
spond to the scattering amplitude and the ra- 
pidity, respectively. In Fig. 1, we show the so- 
lutions to eq. (4) with different initial condi- 
tions at t — 0, together with the correspond- 
ing solutions to the linear equation (3). We 
can learn much from this result. First of all, 
at early time t ~ 0, the solution shows expo- 
nential growth as in the linear case. However, 
as N(t) grows, the nonlinear term (~ N 2 ) 
becomes equally important, and the speed of 
growth is reduced. Eventually at late time, 
the solution approaches to a constant which 
is determined by the asymptotic condition 
dN/dt = 0. This corresponds to the satura- 
tion and unitarization of the gluon number. 
Next, note that two solutions of the logis- 
tic equation with different initial conditions 
approach to each other, and converge to the 
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Figure 1. Solutions to eq. (4) [solid lines] and eq. (3) 
[dashed lines], with different initial conditions. 

same value, while deviation of two solutions 
of the linear equation expands as time goes. 
Namely, in the logistic equation, the initial 
condition dependence disappears as t — > oo, 
which is related to the universality of the BK 
equation. This analogy may look oversimpli- 
fied, but actually most of the phenomena seen 
in the BK equation can be identified in this 
simple example, as far as energy dependence 
of the solution is concerned. 

3 Saturation scale 

Let us turn to the transverse dynamics which 
was ignored above. The most important 
quantity associated with the transverse dy- 
namics is the saturation scale Q s , which gives 
the border line between the saturated (non- 
linear) and non-saturated (linear) regimes. 
Physically, it corresponds to (inverse of) the 
typical transverse size of gluons when the 
transverse plane of a hadron or nucleus is 
filled with gluons. More precisely, it is given 
by the scale Q 2 satisfying 

xG{x,Q 2 ) 



a • p = 1, with a = 



Q 2 ' P nR 2 
Here a is the cross section between the ex- 
ternal probe and a gluon, and p is the gluon 
density per unit transverse area of the tar- 
get with its radius R. Note that the satura- 
tion scale depends upon both a; and atomic 
number A. Its x and A dependences can 
be precisely determined by the BK equation 
or the linear BFKL equation with saturation 



(5) 



"Indeed, the logistic equation (4) is obtained in the 
homogeneous approximation of the FKPP equation, 
which is essentially equivalent to the BK equation. 4 
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boundary condition. When x is not too small, 
one finds 5 ' 6,7 ' 8 

Ql(x, A) oc A 1/3 x~ x (6) 

where A is given by 4.88a s from the LO- 
BFKL, and A - 0.3 for x = 10~ 2 - 10~ 4 from 
(resummed) NLO-BFKL. Also, for the run- 
ning coupling case, A dependence of Q s dis- 
appears at very large energy (small x). 9 This 
particular dependence upon x and A leads 
to an interesting observation that the satu- 
ration scales for the deep inelastic scattering 
(DIS) at HERA and for the Au-Au collisions 
at RHIC are of the same order. Therefore, 
if one finds saturation effects in the HERA 
data, then there is enough reason to expect 
similar things in the RHIC data. 

4 Geometric scaling 

Geometric scaling 10 is one of the most signif- 
icant experimental facts which provides an 
indirect evidence of the CGC. This is a new 
scaling phenomenon at small x meaning that 
the total 7*-proton cross section a^ ^ al in DIS 
depends upon Q 2 and x only via their spe- 
cific combination £ = Q 2 R 2 (x), withi? 2 (x) oc 
x A , A ~ 0.3. Recall that 1/Q is the trans- 
verse size of gluons measured by the virtual 
photon, and, particularly, 1/Q S is the gluon 
size when the transverse plane of the proton 
is filled with gluons. Therefore, if we identify 
Rq(x) = 1/Q s (x), then the scaling variable 
turns out to be a ratio between the trans- 
verse areas of these gluons = o- g /a sat . 
Moreover, if gluons are distributed homoge- 
neously over the transverse plane, which is 
realized in the saturated regime, £ _1 corre- 
sponds to the number of covering times. For 
example, when £ = 1, the transverse plane is 
covered by gluons once which, however, can 
be realized by various size of gluons. Since 
what matters in the saturated regime is the 
effective number of overlapping, it is natural 
that different kinematics can give the same 
value of cross section if the number of cover- 
ing is the same. This is nothing but the geo- 
metric scaling! Furthermore, it is remarkable 
that the theoretical calculation for the x de- 
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Figure 2. "Phase diagram" of a proton in DIS. 

pendence of Q s (x) is consistent with the ex- 
perimentally determined one Q 2 (x) oc a; -0 ' 3 . 

The scenario discussed so far holds natu- 
rally in the saturation regime. On the other 
hand, the data show the geometric scaling up 
to higher values of Q 2 (~ 100 GeV 2 ), which 
is the non-saturated regime described by the 
linear BFKL equation. Indeed, this is again 
well understood in the context of CGC. The 
geometric scaling appears if one solves the 
BFKL equation with a saturation boundary 
condition. As we depart from the saturation 
line toward the linear regime, the effect of 
saturation becomes weaker and weaker and 
eventually disappears. In fact, one can deter- 
mine the window for the geometric scaling: 6 

Ql(x) < Q 2 < Qt(x)/A 2 QCD . (7) 

Hence we recognized that there is a qualita- 
tively different regime in between the CGC 
and DGLAP regimes, as is shown in Fig. 2. 

5 The CGC confronts experiments 

Our understanding of the CGC has been 
developed with the experimental results at 
HERA and RHIC, as we already saw above. 
Below I briefly explain some of the attempts 
to describe/understand the other experimen- 
tal data from the viewpoint of CGC. 

5.1 HERA small x data 
DIS at small x is the cleanest process for mea- 
suring the saturation effects in gluon distribu- 
tion of the proton. Starting from the pioneer- 
ing work by Golec-Biernat and Wusthoff, 12 
there are several efforts to describe the HERA 
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DIS data such as the F 2 structure function 
at small x, in the context of gluon satura- 
tion. Here let us explain the "CGC fit" 13 
which is one of the most successful fits of the 
small x data based on QCD. The CGC fit is 
constructed so as to contain two approximate 
solutions to the BK equation, which are valid 
in the saturation and BFKL regimes, respec- 
tively. In particular, the solution in the lin- 
ear BFKL regime shows the geometric scaling 
and its (small) violation. With only three pa- 
rameters, the CGC fit provides a very nice fit 
for the F 2 structure function with x < 10~ 2 
and in 0.045 < Q 2 < 45 GeV 2 . Meanwhile, it 
turned out that this fit works reasonably well 
even for other observables such as F$ lff and 
the vector meson production. 

5.2 RHIC Au-Au data 

The CGC provides the initial condition for 
the heavy ion collision. Information of the 
initial state could still be seen in the final ob- 
served data. It should be noticed that most of 
the produced particles have small momenta 
less than 1 GeV which is of the same order as 
Q s in RHIC. This observation suggests that 
effects of saturation may be visible in bulk 
quantities such as the multiplicity. Indeed, 
the CGC results 11 for the pseudo-rapidity 
and centrality dependences of the multiplic- 
ity are in good agreement with the data. 

5.3 RHIC d-Au data 

The CGC has recently come under the spot- 
light with the experimental results 14 for the 
nuclear modification factor in the deuteron- 
Au collisions measured by the Brahms exper- 
iment at RHIC. The data show enhancement 
of the ratio at mid-rapidity (the Cronin ef- 
fect) and suppression at forward rapidities. 
In fact, the global behavior of the data is 
qualitatively consistent with the predictions 
made by the CGC. 1546 More recently, de- 
tailed analysis of the ratio was done 17 and it 
has been clarified that the Cronin effect is due 
to the multiple collision and re-distribution 
of the gluons which is properly described by 



the McLerran-Venugopalan model, and that 
the high p± suppression is induced by the 
mismatch of the evolution speed between the 
proton (deuteron) and the nucleus. The nu- 
cleus is closer to saturation and evolves slower 
than the proton. 
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